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Abstract. - In this short note, we present supporting evidence for the replica symmetric 
approach to the random bond g-state Potts models. The evidence is statistically strong enough 
to reject the applicability of the Parisi replica symmetry breaking scheme to this class of models. 
The test we use is a generalization of one formerly proposed by Dotsenko et al. jij and consists 
in measuring scaling laws of disordered-averaged moments of the spin-spin correlation functions. 
Numerical results, obtained via Monte Carlo simulations for several values of q, are shown to 
be in fair agreement with the replica symmetric values computed by using perturbative CFT 
|l| Q for the second and third moments of the q = 3 model. RSB effects, which should increase 
in strength with moment, are unobserved. 



Since its first application to the study of glassy systems, the replica approach has been a 
useful tool in the disordered models research field. However, as was rapidly observed, its 
straightforward application (that is, assuming all replicas are identical) to physical systems 
can lead to serious problems, the most notorious being the negative entropy it gives for the spin 
glass model. These non-physical conclusions can be avoided by breaking the replica symmetry. 
The way in which the symmetry has to be broken in the spin glass problem was first explained 
by Parisi || . However, there are systems in which the symmetry is not broken and where the 
replica symmetric (RS) approach is valid. There is no straightforward arguments to decide 
whether a system studied using the replica method exhibits replica symmetry breaking (RSB) 
or not. 

The aim of this short note is to implement a test that can reveal the presence of RSB 
in disordered local-interaction spin systems. We shall consider disordered Potts models, 
where disorder is introduced via randomness in bond strengths. This problem was studied 
perturbatively using the replica technique j|, 0. There, critical exponents were computed 
assuming RS and the obtained values were shown to be in agreement with Monte Carlo data. 
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However, a possible RSB would not affect those quantities significantly and the apparent 
agreement couldn't rule out its existence. 

Recently, Dotsenko et al. ]l| proposed a test of RSB for random bond Potts model. There, 
they showed that eventual RSB effects could be observed if one considered the disorder- 
averaged moments of spin-spin correlation functions. They studied the second moment and 
agreement was found with RS solutions, thus rejecting RSB-induced deviations that would 
have been greater than the obtained statistical resolution. 

In a recent letter, we exposed perturbative CFT computations of the p-th moment in 
the replica symmetric case. This naturally proposed the search for signs of RSB in these 
correlators. Although we have not computed explicitely the deviations (this would involve a 
tensorial formulation of the Parisi scheme), one can easily convince himself that they should 
become more and more important as p increases. Performing Monte-Carlo simulations, we 
computed the third moment for 3, 4 and 8-state Potts Model. For all these models, we observe 
scaling laws, thus showing that there is no RSB. The associated critical exponents are shown to 
be in fair agreement, for the 3-states model, with CFT predictions, although the perturbative 
expansion is not expected to be very precise for q = 3. This is the first validation of our formula, 
previously exposed [0 , which goes a order further in perturbation than the one originally given 
by Ludwig g. 

Perturbative CFT results. - We shall not repeat here the renormalization group computa- 
tions of higher moments, which can be found, although not in details, in references Jl], ||. We 
rather give a short overview, only stating relevant results. 

The partition function of the nearly-critical (/-states random bond Potts model, is well 
known to be of the form 

Z{f3) = Tr exp{-# - JT X }, (1) 

where Ho is the Hamiltonian of the conformal field theory corresponding to the q-states Potts 
model with coupling constant Jo the same for each bond. The Hamiltonian Hi, being the 
deviation from the critical point induced by disorder is of the form 

Hi = J d 2 xr(x)e(x), (2) 

where t{x) ~ j3J{x) — /3 c Jo is the random temperature parameter. The theory is defined on 
the whole plane. We shall assume, for simplicity, that t(x) has a gaussian distribution for 
each x, with 



t(x) = t = — — (3) 

Pc 

(t(x) ~ t )(t(x<) - t q ) = g 8^(x-x'), (4) 

The usual way of averaging over disorder is to introduce replicas, that is, n identical copies 
of the same model, for which: 

{n r. n \ 

-Y, H o a) - d 2 xT(x)Y,e a (x) • (5) 
a=l J a=l ) 

Taking the average over disorder by performing gaussian integration, one gets 



!n » n » n 

-Y, H o a) - T ° / d 2 xY,£a(x)+g Q / d 2 xY,£a(x)e b (x) } . «i) 
a=l J o=l ^ a^b 
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This is a field theory of n coupled models with coupling action given by 

/n 
d 2 x^e a {x)e b {x). (7) 

a=£b 

Only non-diagonal terms are kept since diagonal ones can be included in the Hamiltonian Hq. 
Moreover, they can be shown to have irrelevant contributions, since their OPE consist of the 
identity plus terms that are irrelevant at the pure fixed point. We now turn our attention to 
the p-th moment of the spin-spin correlation function (a(0)a(R)) p . In terms of replicas, it can 
be written as 

(a(0)a(R))P= lim^^:/ £ a ai (0) • • • a ap (0) £ a bl (R) • • • a bp (R)\ (8) 

Thus, the operator to be renormalized is 

O p (x) = a ai (x)a a2 (x) ■ ■ ■ a ap (x) (9) 

oi 7^ a-2 ■ ■ ■ 7^ CLp, 1 _ o-i _ n, (10) 

perturbed by the interaction term 

6 p (x) = O p exp{-H int } (11) 

= 0p(l-^int + ^(^int) 2 -•••)■ (12) 

Renormalization group computations lead to the identification of a non-trivial fixed point, at 
which we are able to compute the correlation functions. Using scaling laws, we get 

ai=ia 2 ---^ia p 

~ R 2 P A„ ■ 'y 16 ' 

The final result is obtained by using the fixed point value Z(£ R ) ~ e 7 *« R = R 1 * . The 
RG study introduces a parameter e, which can be seen as proportional to the central charge 
deviation of the pure model from the Ising value of 1/2. For the 3-state Potts model, e = 2/15. 
For generic e, one gets (in ||, a should be replaced by —a): 

1 



(^MR)) p ~ (i4) 

with 

A' CTP =A ff - 7 »(p), (15) 



32 



12 



24 



(16) 



and 



J ftT = 61og2 a = 33 



29V3tt 



(17) 
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Thus, perturbed conformal field theory predicts, for the 3-state Potts models, the following 
values for the second and third moments: 

2A; 2 = — - 0.0314 = 0.235 (18) 
1 5 

2A^3 = — - 0.0466 = 0.220 (19) 
1 o 

Monte Carlo Simulations. - To search for signs of RSB, and, in the absence of it, to 
confirm RS values, we performed Monte Carlo simulations of the random bond g-Potts model 
for q — 2,4,8. The method used follows the one in Q. To study scaling effects on the 
correlation functions, we studied square lattices of side L ranging from 10 to 500. Since we 
wanted to exhibit a possible break of the replica symmetry, the algorithm has to be chosen 
in such a way that it doesn't assume the symmetry a priori. For this reason, we simulated 
three configurations of the g-Potts model with same disorder, but different initial conditions 
and independent thermalizations. We computed the products of magnetization 

i=l, L 2 

and 

Q* = J2 <0<^>, (21) 

i=l, L 2 

with (erf) being the thermal average of the local magnetization 

< = 3? ■ ™\ (22) 
where m a is the mean magnetization of lattice a; 

™ a = h E (23) 

i=l,L 2 

It is rather obvious, since all lattices are thcrmalizcd independently, that Q3 and Qi are indeed 
the same as 



^ - t? E (°?°ta> (2 4 ) 

i=l,L 2 

Q* = J2 E (°?°i>- (2 5 ) 



-l,L 2 

Measurement were performed on square lattices with toroidal boundary conditions. The 
Hamiltonian of the simulated model is 

H = - J2 J v + + 5 °h°?) . (26) 

where we took the coupling between nearest neighbours to be 



Jij — Jo or Ji 



(27) 
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with equal probabilities. This makes it possible to make the system self-dual by tuning the 
temperature so that the relation 

= e~^ (28) 



1 + (q - l)e-P J ° 

is obeyed. We chose Jq/J\ — 10 for the simulations with q = 3,4, which is strong enough to 
avoid cross-over effects For the q = 8 model,we rather chose Jq/Ji — 8.5, again because 

this seems the appropriate value to avoid cross-over and minimize the spread of our data set. 

Autocorrelation times were coarsely evaluated and the statistics ajusted in such a way that 
thermal fluctuations can be ignored (typically, for a single disorder configuration, thermaliza- 
tion period was at least 70 auto-correlation times long and at least 200 measures were taken 
(one every auto-correlation time). To average over disorder, we made measurements for 20 
000 disorder configurations (10000 for q = 8). Doing so, one can extract critical exponents 
straightforwardly : 

Q;=KL-p a '^, (29) 

where if is a non-universal constant. The exponent can then be obtained by taking logarithms. 
The results of our simulations, shown in Figures 1, 2 and 3, clearly support the RS scenario. 

2 /p 

In these figures, we present log-log plots of Q p versus L (p = 2, 3), for the three, four and 
eight-state Potts models. By taking the slopes of these graphs, we can extract 2A^. P (q) (which 
is minus the slope). None of the models presented show significant deviations from scaling 
which should arise if the replica symmetry was broken. 

For the 3-state Potts model, the critical exponents associated to the scaling behaviour are in 
fair agreement with the values predicted by perturbative CFT computations. The deviations 
from the pure model behaviour are: 

2 7 *(2) = 0.0387 Monte Carlo (30) 
= 0.0314 CFT prediction (31) 

2 7 *(3) = 0.0648 Monte Carlo (32) 
= 0.0466 CFT prediction (33) 

The numerical agreement is indeed quite surprising, especially for the third moment, where 
the perturbative expansion is near the end of its validity region Eqj. Olson and Young |8| also 
computed spin-spin correlation functions moments, but in a different optic and with a different 
method. Our values for the exponents, presented in Table I, are in fair agreement with theirs, 
altough they seem to be systematically lower. Using other methods, Palagyi, Chatelain et 
al. [|9j obtain values that confirm this discrepancy. Our values are equal to theirs, within 
statistical errors. 



Conclusion. We believe that the presented evidence is enough to rule out the RSB 
scenario in random bond Potts models. If this symmetry was broken following Parisi's scheme, 
deviations from the observed scaling laws would be, for the second moment, of the order of 
10%, and thus would be easily observed. One can convince himself that the deviation should 
become more apparent for the third moment, something which is clearly not observed. 

It will be interesting to see how more precise numerical methods, such as transfer matrices 
iterations flio[ , could give accurate values for moments via cumulant expansions (for integer 
and non-integer values of q). 
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Fig. 1. - Log- Log plot of Q p with p = 2 (lower line) and p = 3 (higher line) for the random 3-state 
Potts model 
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Fig. 2. - Log-Log plot of Q p P with p = 2 (lower line) and p — 3 (higher line) for the random 4-state 
Potts model 



We would like to thank VI. S. Dotscnko, M. Picco, J.L. Jacobsen, C. Chatelain and B. 
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Fig. 3. - Log-Log plot of Q v with p = 2 (lower line) and p = 3 (higher line) for the random 8-state 
Potts model 
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2A; 2 


2A^ 3 


3 


0.228(1) 


0.202(2) 


4 


0.231(2) 


0.198(2) 


8 


0.229(2) 


0.184(3) 



Table I. - The values of the spin exponents 2A^. 2 and 2A^. 3 , for q = 3,4, 8. The normalisation was 
chosen so that exponents for the same value of q could be compared. The number in parentheses is the 
statistical error in the last decimal place 
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